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Motivated by the recent optical conductivity experiments on Bi2Sr2CaCu2 0s+5 films, we examine 
the possible origin of low-frequency dissipation in the superconducting state. In the presence of 
spatial inhomogeneity of the local phase stiffness p 3 , it is shown that some spectral weight is removed 
from uj = to finite frequencies and contribute to dissipation. A case where both p s and the local 
normal fluid density are inhomogeneous is also considered. We find an enhanced dissipation at low 
frequency if the two variations are anti-correlated. 

PACS numbers:74.20.-z,74.20.-g,74.20.De 



Identifying the nature of quasiparticles and their dy- 
namics is an important part of the overall understanding 
of the properties of matter. In d-wave cuprate super- 
conductors, gap-vanishing nodes are responsible for an 
appreciable density of quasiparticles even for a supercon- 
ducting state. Optical conductivity measurements car- 
ried out on YBCO crystal Q shows that the quasiparticle 
dynamics in the superconducting state obeys the Drude 
form, 



Pn{T)T n {T) 

l + w 2 7*(T)' 



(1) 



where p n (T) is the T-linear normal fluid fraction, and 
r n (T) is the quasiparticle lifetime at a given tempera- 
ture. As the temperature is lowered er„(u>,T) is seen to 
decrease due to vanishing quasiparticle density. 

For Bi 2 Sr 2 CaCu20 8+ 5 (Bi:2212), on the other hand, 
a similar microwave measurement shows substantial 
amount of residual conductivity as u> — > and [§, 
even though the operating temperature and frequency 
ranges of the measurement were comparable to those in 
Rcf . . It was recently shown that a good fit to the op- 
tical conductivity data for slightly underdoped Bi:2212 
thin film (T C =85K) can be achieved if one introduces 
another term, denoted o~ c , in addition to the two- fluid 
contribution in the conductivity 



a c (w,T) = k 
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The low-frequency conductivity o~ c now grows as the su- 
perfluid density p s {T) upon lowering the temperature. 
Further studies showed that k, which is a measure of the 
amount of spectral weight removed from the condensate, 
grows monotonically with increasing doping concentra- 
tion. 

Such behavior of the optical conductivity is unusual 
from the simple two-fluid point of view, and deserves the- 
oretical attention It is also essential to clarify the 



source of discrepancy of the optical data among differ- 
ent materials. In the context of Josephson-coupled su- 
perconducting grains, Barabasi and Stroud showed that 
an inhomogeneous distribution of the Josephson coupling 
strength J in superconducting arrays gives rise to a shift 
of spectral weight under ui — (the condensate frac- 
tion) to finite frequencies with an amount proportional 
to ((5J) 2 )/J, where J is the mean phase stiffness of the 
system [Q. Roughly speaking, the inhomogeneity in J 
gives rise to scattering of superfluid density fluctuations 
(spin waves) which tends to move some of the states from 
the condensate into the dissipative part. 

Evidence in support of inhomogeneous nature of the 
superconducting state in Bi:2212 has been mounting 
from STM images of the BSCCO surface [§-§. Theo- 
retical consideration shows that the randomness of the 
dopant oxygen, inherent in the doping process for Bis- 
muth compounds, does produce an inhomogeneous hole 
distribution in the copper-oxygen plane ||. The zero- 
temperature phase stiffness is in turn determined by the 
local hole density in cuprate superconductors, and is very 
likely to be disordered. The purpose of this paper is to 
calculate the optical conductivity induced by the ran- 
domness in the local phase stiffnes /5 s ( r )- Inherent in 
our model is the assumption that the inhomogeneity is 
in fact present in the bulk, although the direct evidence 
from STM is confined to the surface layer of the com- 
pound. Other evidence in support of the bulk-disordered 
nature of the Bi:2212 superconductors can be found, for 
example, in transport, neutron scattering, and ARPES 
measurements of Bi:2212 against a similar measurement 
for YBCO §. 

Our approach is based on the effective action for the 
phase of the condensate order parameter, <fi. Starting 
from a Lagrangian which gives two-fluid description of 
the low-energy excitations in charged superfluids, we ar- 
rive at an effective action for (f> alone through integrat- 
ing out other degrees of freedom. A general formula for 
the optical conductivity in the presence of phase stiffness 
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disorder is derived via the replica method, and using the disorder strength, 
effective action for <f>, computed to lowest order in the 



Effective action of the phase- fluctuating (2+l)-dimensional superfluid can be written down in the following general 
form 

S =\J d 3 xd 3 y[K(x - y) + S Ps (r)S(x-y)](grsid<j)+A){x){grsid(t>+AXy) 

cPxtfyidot+AoHxp-Hx - y)(d o 0+A Q )(y). (3) 



We denote the space-time coordinates as x and y. All the 
linear-response properties are summarized by the non- 
local kernels K(x — y) and U~ 1 (x — y), which will be 
specified later. As the above action shows, the stiffness 
disorder Sp s is quenched, and does not fluctuate in time. 
We further assume that 5p s (r) obeys the Gaussian dis- 
tribution 



P[Sp s (r)} cx exp 
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The disorder average of the partition function Z can be 
carried out in the usual manner by first re-writing In Z 
as (Z N — 1)/N, then integrating over 6p s (r). This leads 
to a replica action which is N copies of the action in Eq. 
(0) without Sp s , plus the disorder-induced term 
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where m,n are the replica indices. 

The current-current correlation function IT 
from 



(5) 



follows 



Ily (ri-r 2 ,Ti-T 2 ) = - lim 



d 2 {lnZ) 



n^o dAi(riTi)dAj(r 2 T2) ' 
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using the replica action defined above. Also in an optical 
experiment we are interested in the transverse piece of 



EL,-. Contributions of Eq. 



to Hij is two-fold. First, 



expanding Eq. (|^) one obtains (dropping replica indices) 
~ (grad</>) 2 (grad</>) 2 term which, to first order in g, can 
be replaced by ~ ((grad</>) 2 )(grad</>) 2 and modifies the 
kernel K(k). Secondly there is ~ (A • grad^) 2 , which 
will contribute ~ (grad</) grad</>) in the current-current 
correlation function. As a result the following form is 
obtained for the transverse part: 

U t (k) = K{k)-igG(iuj m ). (7) 

is the Fourier transform of the 



Here, K(k), k — (k,iui. 
kernel K(x — y), and 



G(iu> m ) 



d 2 k 
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k 2 (</ ) (-fc)0(fc)). 
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The momentum integral is restricted to |k| < A, which 
is the inverse of the typical length scale of the inhomo- 
geneity. The boson-boson correlator is given by 

(cb(-k)cb(k))- 1 = s(k, iu m ) = k 2 K(k) + uj 2 m /U{k). (9) 

The optical conductivity is in turn obtained from IT t 
through 



Il t (k = 0,iw m ) 



(10) 



Equations (^)-(|T^) form the basis for calculating the 
optical conductivity of a given superfluid in the presence 
of phase stiffness disorder. The properties of a super- 
fluid are encoded in the linear response functions K (k) 
and U~ l (k), which in turn can be derived from a specific 
microscopic model. Once they are determined, optical 
conductivity immediately follows from the above formu- 
lae. 

We model the low-energy dynamics of a superconduc- 
tor as charged two-fluid system described by the action 
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S n + ^(gr a dA ) 2 . 
on 
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In the case of d-wave superconductors, the quasiparti- 
cle action S n is that of Dirac particles. The super- and 
normal-fluid three-currents are respectively denoted as 
Jsp, = (ps—p, Js), and Jnp, (p = 0, 1, 2). One recovers the 
current conservation law, d^(J sfJl + J„ M ) = 0, upon inte- 
grating out the phase 4>. There is no separate conserva- 
tion of J sp or J njl alone since conversion of superfluid into 
normal fluid and vice versa can take place through scat- 
tering processes. The gauge field action e (gradA ) 2 /87r 
generates the long-range Coulomb interaction between 
charges. The bare stiffness and compressibility of the 
superfluid are given by K and w -1 , respectively. 

On integrating out the normal current J nfjL and the 
supercurrent J s ^ successively we arrive at an effective 
action: 
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S=y(grad0+A) 2 + i-(d o < 



-A,) 2 



+ hd^ + A^n^id^ + A^ + ^igradAo) 2 . (12) 

The result of integration over the normal current J nfi 
has been expanded up to quadratic order in d^cj) + A^, 
therefore k^u is the correlation function of the normal 
3-current, {JntiJnv)- We further write koo = Xn, an d 
Kn = ^22 = Kn- For the off-diagonal terms, xoi = 
due to the absence of time-reversal-symmetry breaking, 
and Xi^j = assuming a typical nonmagnetic impurity 
scattering of quasiparticles. 

Let us first review the case of homogeneous supercon- 
ductor. We set K — p s , and write n n in terms of the 
quasiparticle conductivity a n as 



-Pn(T) 



(13) 



With this substitution, and in the absence of disorder, 
the transverse current-current correlation function II t be- 
comes p s — p n (T) + Ionian. Once p s — p n (T) is identi- 
fied as the finite-temperature superfluid density p s (T), 
the conductivity at long wavelength (k = 0) is obtained 
fromEq. @ 



<r(w,T) = 



i-Ps{T) 
u) + id 



(14) 



This is nothing but the two-fluid result for the conductiv- 
ity. We have thus demonstrated that our starting point in 
Eq. ( p"l| ) is indeed consistent with the two-fluid picture. 

Finally upon integrating out the scalar potential Aq in 
action (13), we obtain the boson propagator 



s(k, iu> m ) = {p s (T) + \uj m \a n )l 



1 + 4tt(u- 1 +Xn)Aok 2 ' 



(15) 



Knowledge of this quantity is not necessary for the eval- 
uation of the conductivity in the homogeneous supercon- 
ductor, but will be important for the disorder-induced 
contribution as one can see from Eqs. (@)-(|). One can 
in fact identify p s (T) + |w m |cr ra and u _1 + x„ as current- 
current and density-density correlation functions of the 
combined superfluid + normal system. Inclusion of the 
Coulomb interaction renormalizes the latter quantity in 
the RPA manner [jl0|. We can now proceed using this 
boson propagator to understand the effects of inhomo- 
geneity. To clarify the physical content of the expression, 
Eq. (|l5|), we introduce 



7 = 71 + ij2 = 



47T 



Xn 



and write (after iu m — > u> + i5) 



s(k,b)) = — [^{p s (T)-iwa n )-{uj+i8Y 



(16) 



(17) 



Without <x„, the terms inside the bracket gives the 
Josephson plasma dispersion of the charged superfluid. 
With the quasiparticles present, damping occurs in the 
plasma mode with a characteristic lifetime ~ a n /p s (T). 

When s(k,(j) is substituted in Eq. (0) the disorder- 
induced conductivity becomes 



a' 1 {u) = 



gA 2 
Ps{T) 

-— Im 



6(d) + a[(b)) 



d 2 k 



(27r) 2 u) 2 -"f(p s (T)^iu)a n ) 



(18) 



For convenience we have written down only the real part. 
The amount of spectral weight shift away from the con- 
densate is indeed proportional to (8p 2 s ) / p s (T) since 
g A 2 ~ (8p 2 s ). They will re- appear at non-zero frequen- 
cies in a'^ui). We examine the behavior of cr^(u)) under 
a number of circumstances. 



I) Zero quasiparticle dissipation: Take a n = 0. In this 
case, 



a[ (b)) w -Agu 



d 2 k 
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(27r) 2 [u) 2 - 7lPs (T)} 2 + [ 1 2Ps(T)} 2 - 



(19) 



For small 72, the integrand is narrowly centered around 
the plasma frequency, to 2 ~ liPs, and zero otherwise. 
Hence, the dissipation arises only if u) exceeds a charac- 
teristic value comparable to Josephson plasma frequency. 
The experimental observation of low-frequency dissipa- 
tion H is therefore impossible to explain in this case. 

II) Finite quasiparticle dissipation: We are after all 
interested in frequency ranges much smaller than the 
plasma frequency. In this case, the integrand appearing 
in Eq. (fl8|) is simplified to (iuja n — p s (T)) _1 and 



a[(b)) 



gA 2 



Ps(Tf 



(20) 



This is a Lorentzian centered at uj = 0, with the width 
governed by p s (T)/a n . By integrating this expression 
over positive uj, one finds that all of the spectral weight 
removed from the condensate indeed re-appears in a[(b)). 

We have thus demonstrated the existence of a new 
low-energy dissipation term, Eq. (p0|), for the superfluid 
which was absent in the conventional two-fluid picture. 
Presence of this effect rests on 1) residual quasiparticle 
dissipation (er n ^ 0) and 2) inhomogeneity in the super- 
fluid phase stiffness (gA 2 > 0). Given that the doping 
process almost always results in an inhomogeneous distri- 
bution of holes, we expect our result to hold for a variety 
of other cuprate superconductors. 
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Up to this point it is assumed that a n remains homoge- 
neous throughout the system, while the local superfluid 
density undergoes a spatial fluctuation. It is certainly 
possible, however, that the local normal fluid density is 
disordered as well as p s . For a sufficiently smooth vari- 
ation, we can treat the local value of the quasiparticlc 
conductivity <r„(r) to be proportional to p n (r). In par- 
ticular we will consider a situation where the variation of 
normal fluid density is tied to that of the superfluid by 
Sp n (r) = XSp s (r), with A taking on positive (correlated) 
or negative (anti-correlated) values. The quasiparticle 
conductivity is accordingly 



,(r) 



XSp s (r)T n 



(21) 



We assume that the quasiparticle conductivity remains 
in the d.c. limit, LOT n <C 1. 

Replica averaging over Sp s (r) can be carried out 
straightforwardly as before and yields the following com- 
plex conductivity due to inhomogeneity: 



—i gA 2 (1 — i\T n Lu) 2 
u)+iS 7r p s (T) — ia n u> 



(22) 



Our previous discussion belongs to A 
is now given by (u> > 0) 



0. The real part 



cti (u>) 



gA 2 a n - 2X Ps {T) Tn 
tt p s (T) 2 + a 2 n u 2 ' 



(23) 



ignoring a {ujT n ) 2 correction. Compared to our previ- 
ous result in Eq. (pp[), one finds an increase in the 
low- frequency conductivity in the case of anti-correlated 
spatial variations of the normal and superfluid densities 
(A < 0), and a decreased conductivity otherwise. A sim- 
ilar conclusion has been reached in a recent preprint by 
Orenstein ||. 

We now discuss how the formulas derived above, in 
particular Eqs. ( pfj| ) and (|2lj|), compare with the exper- 
iment. The damping rate t^ 1 = p s /a n turns out to be 
100K w 2THz, based on the zero-temperature estimate 
of p s and o~ n from the microwave and terahertz data 
||. This value is consistent with the measured value 
of 1.5THz f|. If we re-cast Eq. @ in the form of 
Eq. (^|) using the definition t c = o~ n /p s (T), we arrive at 
k = gA 2 /irp s (T) 2 ~ (Sp 2 ) / p s (T) 2 . The significance of k 
is thus clarified, as the relative fraction of the inhomo- 
geneity present in the superfluid density. Furthermore, 
re-writing Eq. (j23|) in the same form gives 



k(A) = k(A = 0) x 1 



2Ap s (T) 



Pn 



(24) 



The T — > estimate of the superfluid and normal fluid 
densities turn out to be comparable in available sam- 
ples H , and the degree of correlation A between normal 



fluid and superfluid densities turn out to play an impor- 
tant role in estimating n. Experiments carried out on 
Bi:2212 films of varying doping concentrations ranging 
from underdoped (T C =51K) to overdoped (T C =75K) in- 
dicate that k monotonically increases from 0.08 to 0.46 
H with doping. Equation ( p4| ) suggests that one way to 
understand the observed behavior of K is to say that the 
superfluid and normal fluid densities in the ground state 
become more and more anti-correlated upon over-doping. 

In summary, we have considered the optical conduc- 
tivity response of a charged superfluid in the presence 
of inhomogeneous phase stiffness distribution. In the ho- 
mogeneous system only the quasiparticle excitations con- 
tribute to low-frequency dissipation. On the other hand, 
we find that some spectral weight is removed from the 
condensate and shows up as low-frequency dissipation 
due to spatial fluctuations of the local phase stiffness. In 
the model we considered, the low-frequency distribution 
follows a Lorentzian. The spectral weight lost at to = is 
proportional to {8p 2 s ) / p s {T). In the case where the nor- 
mal fluid density p n is also disordered, the low-frequency 
dissipation is still present, and is increased relative to the 
Sp n = case provided Sp n (r) oc — Sp s (r), i.e. they are 
anti-correlated. A decreased conductivity is found if they 
are spatially correlated. 
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